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ABSTRACT 

We extend the Halo Occupation Distribution (HOD) framework to generate mock galaxy 
catalogs exhibiting varying levels of “galactic conformity”, which has emerged as a poten¬ 
tially powerful probe of environmental effects in galaxy evolution. Our model correlates 
galaxy colours in a group with the concentration of the common parent dark halo through 
a “group quenching efficiency” p which makes older, more concentrated halos at fixed 
mass preferentially host redder galaxies. We find that, for a specific value of p, this 1-halo 
conformity matches corresponding measurements in a group catalog based on the Sloan 
Digital Sky Survey. Our mocks also display conformity at large separations from isolated 
objects, potentially an imprint of halo assembly bias. A detailed study - using mocks with 
assembly bias erased while keeping 1-halo conformity intact - reveals a rather nuanced 
situation, however. At separations < 4Mpc, conformity is mainly a 1-halo effect domi¬ 
nated by the largest halos and is not a robust indicator of assembly bias. Only at very 
large separations (> 8Mpc) does genuine 2-halo conformity, driven by the assembly bias 
of small halos, manifest distinctly. We explain all these trends in standard Halo Model 
terms. Our model opens the door to parametrized HOD analyses that self-consistently 
account for galactic conformity at all scales. 

Key words: galaxies: formation - cosmology: dark matter, large-scale structure of Uni¬ 
verse - methods: N-body, numerical 


1 INTRODUCTION 


The processes that cause star formation activity in galaxies 
to diminish or “quench” are of great interest in understand¬ 
ing galaxy evolution. A number of observational and (semi- 


)numerical studies (van den Bosch et al.|2008||Ross & Brunner 

200911 Wetzel & White 

2010 

van der Wei et al.|2010 [Hopkins 

et al.|2010|| 

More et al. 

2011 

Prescott et al.|2011t Wetzel et al. 

2 OI 2 I 2 OI 3 I 

Gen|2014 

lirschmann et al.|2014 1 have attempted 


to distinguish between the relative importance of galaxy merg¬ 


ers (Toomre & Toomre 19721 and wider halo-level physical 


processes such as ram-pressure stripping of gas {Gunn & Gott 
19721, “strangulation” ( Larson et al.||1980 Balogh & Morris 
20001, “harassment” ( Moore et al.||1996| ), etc. that may lead 
to quenching in group environments (for a recent review, see 
Somerville fc Dave|2014 1. Even where halo-level effects are in¬ 
dicated, there is ongoing debate over whether the location of 
a galaxy within its group (i.e. whether it is a “central” galaxy 
or a “satellite”) is paramount ( Baugh|[2006 Peng et al.||2012 
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Kovac et al. 20141, or whether the nature of the group as a 


whole is more important in quenching its member galaxies 
(Osmond fc Ponman|200^ van den Bosch et al.|2008 Carollo 


et al.|20T4 Hartley et al.|20T5 Knobel et al.||20l4 l. 


In this context, a particularly interesting phenomenon 
that has received attention recently is that of “galactic con¬ 
formity” ( [Weinmann et al.|2006p . This is the observation that 
satellite galaxies in groups whose central galaxy is quenched 
(or red) are preferentially quenched, even when the groups are 
restricted to reside in dark matter halos of the same mass (see 
[Knobel et al.|20T^ for an in-depth study extending this state¬ 


ment to a variety of measures of environment). Weinmann 


et al.[( 

2006 

et al. 

2005 


20051 constructed using galaxies in the Sloan Digital 


Sky Survey (SDSS). More recently, a similar effect has been 
observed by [Kauffmann et ah] ( |2013| ) when studying the star 
formation rates (SFR) of SDSS galaxies at large spatial sepa¬ 
rations (< 4Mpc) from isolated objects. 

The observation of galactic conformity contradicts an as¬ 
sumption typically made in halo occupation models of galaxy 
distributions, namely, that galaxy properties such as luminos¬ 
ity and colour distributions are determined solely by the mass 
of the parent dark halo of the group to which the galaxy 
belongs. Modelling conformity then becomes important in 
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2 A. Paranjape, et al. 


order to better understand the galaxy-dark matter connec¬ 


tion (Zentner et al. 20141, with consequences for both semi- 


analytical models of galaxy evolution | Wang et al. 2013 Kauff- 
mann et al.||2013l as well as statistical studies aimed at con¬ 


straining cosmological parameters (van den Bosch et al.|20f3 


More et al.|2015 Coupon et al.|20T5l. 

It is tempting to ascribe a physical origin to galactic con¬ 


formity, as Hartley et al. (20151 do. Since the seminal work 
of Fabian et al. ( |2003 1, evidence for regulation of the hot gas 
temperature in clusters of galaxies, by the super-massive black 
hole (SMBH) at the centre of the brightest cluster galaxy, has 
been growing ( Croton et al.|2006 Zhuravleva et al.|2014 1 . Al¬ 
though it is not clear whether a similar regulation exists for 
group-mass halos, a coupling of the star-formation properties 
of different group members via the hot gas content of their halo 
provides a natural explanation for conformity; a coupling that 


could even be established many Gyrs previously (Rawlings & 
Jarvis|2004 1. A coupling such as this may also exist due to re¬ 


processing of the gas content of individual galaxies in a group 
( Oppenheimer et al.|2010 Birrer et al. ]M4| |. However, this is 
far from a unique solution. We might equally speculate that 
the virial shock heating of infalling gas depend on the halo’s 
location within the cosmic web, perhaps being more effective 
in hlaments than in voids at low halo masses for instance. In 
each of these hypothetical cases the conformity signal found by 


Weinmann et al. (20061 and others is driven by a hidden vari¬ 


able (past/present SMBH activity or large scale environment), 
i.e. a quantity that was not accounted for in those studies. 
Clearly identifying the correct hidden variable(s), those that 
“explain” conformity, will lead us to a greater understanding 
of the physical nature of quenching in galaxy groups. 

An obvious suspect as to the hidden variable, and there¬ 
fore the possible origin of galactic conformity, is the assembly 
history of the parent halo. The phenomenon of halo assem¬ 
bly bias, in which low mass older halos tend to cluster more 
strongly at scales > lOMpc than younger, more recently as¬ 
sembled halos of the same mass (while this trend reverses for 


more massive halos) is well established ( 

Sheth & Tormen 

2004 

Gao et al. |2005| 

Wechsler et al.| 

2006 

Jing et al.||2007 

Des- 

jacques|20081 |Da 

al et al.|2008| F 

ahn et al.|20091|Fakhouri & 

Ma[[2010 1. Halo age also correlates well with other halo prop- 


erties - we will focus on halo concentration which correlates 
positively with age (Navarro et al.|lSi^ Wechsler et al.|2002 I 


- and halo assembly bias has been observed to extend to these 
properties as well; e.g., at fixed low mass, more concentrated 
halos typically cluster more strongly than less concentrated 
ones. 

Viewed from the present day, a halo that has a low late¬ 
time mass accretion (with respect to other similar mass haloes) 
would have a higher formation redshift and, statistically, a 
greater concentration. In relatively low-mass haloes the build 
up of dark matter mass and the accretion of baryons by the 
central galaxy are quite tightly coupled, because the time 
scale for cooling to offset gravitational heating is shorter than 


the free-fall time White & Rees (19781. The galaxies hosted 


by early-forming haloes would therefore experience a reduced 
supply of star-forming gas with respect to the larger galaxy 
population, and in the most extreme cases could result in a 
higher probability of being quenched. A plausible hypothesis, 
then, is that halo assembly bias affects galaxy formation and 


leads to a galaxy assembly bias (Tinker et al. 20121, which 


leaves an imprint in the form of galactic conformity (Hearin 
|et al.|[2015[). More specifically, according to this hypothesis 


galactic conformity both within individual groups as well as 
at large scales may be explained if, at fixed halo mass and 
galaxy luminosity/stellar mass, quenched or red galaxies pref¬ 
erentially reside in older (sub)halos. 

One point of interest in this regard is that the conformity 


signal observed by Kauffmann et al. (20131 at scales < 4Mpc is 


quite strong (roughly an order of magnitude difference in the 
median specific SFR of galaxies surrounding isolated objects 
that belong to the upper and lower 25th percentiles of specific 
SFR). Halo assembly bias, on the other hand, has been argued 
to be relatively unimportant in modelling, e.g., the dependence 
of galaxy 2-point correlation functions on properties such as 
luminosity, colour and SFR at similar scales (e.g., [Skibba fc| 
Sheth|[2b09 Deason et al.|[2013|. It is therefore important to 


ask how assembly bias can be simultaneously consistent with 
these observations, or whether some other mechanism could 
generate the conformity signal. 

In principle, such a test can be performed by embedding 
an assumed conformity-inducing correlation between galaxy 
and halo properties in an analytical Halo Model framework 
(for a review, see Cooray fc Sheth|2002 l, since it is feasible to 
accurately model property-dependent clustering in this lan¬ 
guage ( [Wechsler et al.||2006| ). However, comparing analytical 
calculations directly with observations can become quite in¬ 
volved in the presence of possibly complex selection criteria 
used in dehning the observed sample and the signal itself. A 
more efficient method, then, is to use V-body simulations of 
dark matter and construct mock galaxy catalogs in which one 
has full control over switching the conformity signal on and off; 
these mocks can then be analysed in the same way as the ob¬ 
served sample for a fair comparison. Such mock catalogs have 
recently been presented by [Hearin fc Watsonj ( [2013[ ); [Masaki[ 
et al. (2013b I; Hearin et al. (2014a[) who model galactic con¬ 


formity by introducing a rank ordering of galaxy colours or 
SFR by suitably defined measures of (sub)halo age. The “age 
matching” mocks of Hearin fc Watson[( 20131 and Hearin et al. 


(2014a I have been shown to successfully reproduce a number 


of observed trends such as the 2-point correlation function and 


galaxy-galaxy leasing signal of SDSS galaxies (Hearin et al. 


2014bI, the radial profile of the satellite quenched fraction in 
groups ( [Watson et al.[[2015p , etc., in addition to containing a 
conformity signal. The rank ordering in these mocks leads to a 
fixed strength of galactic conformity, however (see, e.g., Figure 
3 in [Hearin et al. ]20T^ . In the absence of conclusive evidence 
as to the nature of the hidden variable that explains confor¬ 
mity, it is then interesting to explore alternative algorithms, 
particularly those that implement a tunable effect which can 
then be compared with observations. 

Another relevant issue is that it is important to be able 
to distinguish cleanly between conformity within individual 
groups and conformity effects due to spatial correlations be¬ 
tween distinct halos (these were respectively dubbed “l-halo” 
and “2-halo conformity” by Hearin et al. 20151. This is because 


the distinction between what is inside and outside a dark halo 
becomes fuzzy in any analysis that averages over a large range 
of halo masses. E.g., we would like to ask the following: out 
to what spatial separations might one see the effects of con¬ 
formity in a model that has only 1-halo conformity built in 
and does not know about the large scale environmental corre¬ 
lations due to halo assembly bias? This might plausibly be the 
case in a halo-specific gas regulation mechanism that couples 
to star formation activity as mentioned above. This question 
is particularly interesting in the context of the apparent con- 
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Conformity from correlating colour and concentration 3 


flict in the expected strength of galaxy assembly bias and the 
observed strength of galactic conformity mentioned above. 

In this paper we introduce an algorithm that allows us 
the level of control needed to address the above issues. Our 
algorithm is a modihcation and extension of the one described 
in Skibba & Sheth (2009 hereafter, S09) to model galaxy posi¬ 


tions, luminosities and colours; we model galactic conformity 
by introducing a positive correlation between galaxy colour 
and the concentration of the parent dark halo of the group to 
which the galaxy belongs (i.e., we work at the level of groups 
and not subhalos). The strength of this correlation is deter¬ 
mined by adjusting the value of a parameter p which we inter¬ 
pret as a “group quenching efficiency”, the terminology being 
motivated by similar quantities studied by earlier authors (see, 
e.g.,|van den Bosch et al.|2008[[Peng et al.|2010||Knobel et al.| 

2014[ ). In order to answer the question regarding large scale 
effects raised above, we also explore a model where halo con¬ 
centrations are randomized among halos of fixed mass before 
correlating them with the galaxy colours, thereby erasing any 
large scale clustering of the colours due to halo assembly bias 
while keeping average group-specific properties (including 1- 
halo conformity) intact. 

The paper is organised as follows. In section[^we describe 
the S09 algorithm, followed by our modifications for correlat¬ 
ing galaxy colours with halo concentrations and calculating 
stellar masses. In section we describe the A^-body simula¬ 
tions on which our mock galaxy catalogs are based. Section]^ 
describes the SDSS-based group catalog that we use for com¬ 
parison with our mocks, together with various fitting functions 
derived from this sample which inform our mocks. Our results 
are described in section followed by a discussion in sec¬ 
tion We conclude in section We will use SDSS galaxy 
properties K-corrected to redshift « = 0.1, usually denoted 
by a superscript, e.g., ^'^r for the r band. Throughout, we 
will denote Mo.i^ — 51ogjo(^) as Mr, ^'^{g — r) as g — r and 
quote stellar masses m* in units of h~^MQ and halo masses 
m in units of where Hq = 100 hkms“^Mpc“^ is the 

Hubble constant. We use a fiat ACDM cosmology with param¬ 
eters fim = 0.25, fib = 0.045, h = 0.7, ug = 0.8 and Ua = 0.96, 
which are consistent with the 5-year results of the WMAP 
experiment (Hinshaw et al.|2009|. 


2 ALGORITHM 

Our basic algorithm is borrowed from S09 and is itself an ex¬ 
tension of the algorithm described by Skibba et al. (20061 to 
include galaxy colours. The algorithm uses a Halo Occupation 
Distribution function (HOD; Berlind & Weinberg 20021 cal¬ 
ibrated on SDSS luminosity-dependent projected clustering 
measurements to create a luminosity-complete mock galaxy 
catalog. We describe this algorithm and our modifications be¬ 
low. 


2.1 The S09 algorithm 

2.1.1 Luminosities, positions and velocities of centrals and 
satellites 

The S09 algorithm explicitly implements the so-called central- 
satellite split. A fraction /cen(< Afr,max|m) of m-halos (i.e., ha¬ 
los with masses in the range (m, m -|- dm)) is chosen to have a 
central galaxy brighter than the luminosity threshold Mr,max- 


Each m-halo with a central is then assigned a number of satel¬ 
lites drawn from a Poisson distribution with mean Naa.t{< 
Mr,max I m). The luminosities of centrals and satellites are 
then assigned using the distributions /cen(< Mr|m)//cen(< 
Mr,max I m) and Naa.t{< Mr | m) /Agat (< Mr,max|m), respec¬ 
tively. The functions /cen(< Mr|m) and Naa,t{< Mr|m) de¬ 
fine the HOD, with the mean number of galaxies brighter 
than Mr residing in m-halos given by (Agai|m) = /cen(< 
Mr|m) [l -I- Aaat(< Mr|m)] . In this work we will use the forms 
calibrated by Zehavi et al.| | |2011 1 ; specifically, we use an in¬ 
terpolation kindly provided by Ramin Skibba. Our fiducial 
cosmology has the same parameter values as used by |Zeha^ 


et al. (20111 for their HOD analysis. Satellite luminosities are 


assigned after the central ones, and the algorithm ensures that 
the central is the brightest galaxy of the group. We will return 
to this point later. 

Each central is placed at the center-of-mass of its par¬ 
ent dark matter halo and is assigned the velocity of the halo. 
The satellites are distributed around the centrals according 


to a truncated Navarro-Frenk-White (NEW) profile (Navarro 


et al. 19961, using a halo concentration as described below. 


and are assigned random velocities relative to the central that 
are drawn from a Maxwell-Boltzmann distribution that scales 
with halo mas^ This procedure for assigning satellite posi¬ 
tions and velocities is more efficient, but perhaps less accurate, 
than identifying satellites with subhalo positions (e.g., there 
is no information regarding infall times and number of peri- 
center passages, but we need not run a very high resolution 
simulation and track low mass subhalos). 


2.1.2 Galaxy colours and the central-satellite split 


Having assigned galaxy positions, velocities and luminosities, 
galaxy colours are assigned by drawing from double-Gaussian 
fits to the observed distribution of g — r colours at fixed lumi¬ 
nosity. The two components of the double-Gaussian - ‘red’ and 
‘blue’ - have means, variances and relative fraction as func¬ 
tions of luminosity that are fit to data. We give the results of 
such fits to SDSS data in section IT^ Since these fits work at 
the level of the full galaxy sample, there is some freedom in 
deciding what fraction of centrals and satellites at fixed lumi¬ 
nosity must be labelled ‘red’ with the above dehnition. This 
is a somewhat subtle issue which merits some discussion. 

In general, one might assume that the central and satel¬ 
lite red fractions depend on both luminosity of the object as 
well as the mass of the parent halo, and we would denote these 
quantities as p(red|cen, Mr, m) and p(red|sat, Mr, m), respec¬ 
tively. S09 showed that a model in which these quantities do 
not depend on halo mass (but have a non-trivial luminosity de¬ 
pendence) is consistent with measurements of colour-marked 
clustering of SDSS galaxies. However, colour-dependent clus¬ 
tering measurements will also depend on the level of confor¬ 
mity, which is what we are trying to model here (and which 
S09 did not do). We therefore adopt a simpler approach: we 
continue to assume that these red fractions are independent 
of halo mass, but fix the form of p(red|sat, M^) by demand¬ 
ing agreement with a direct measurement of this quantity in 


an SDSS-based group catalog (Yang et al. 2007 henceforth. 


^ We note that our analysis in this paper does not use information 
regarding galaxy velocities; this will however be essential in future 
work when comparing, e.g., to observations of clustering in redshift 
space. 
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Y07). The form of p(red|cen, Mr) is then fixed by the assumed 
HOD, while the aZ/-galaxy red fraction p(red|Mr,m ) inh erits 
a mass dependence from the HOD (see Appendix |Al[ P] In 
Appendix | A1 1 we also explore the consequences of a halo mass 
dependence in the red fraction of centrals. 

The S09 algorithm has been used in several studies of 
galaxy environments (e.g., Muldrew et al.|[2012 Skibba et al. 


20131 and the luminosity and color dependence of galaxy clus¬ 


tering (e.g., Skibba et al. 2014 Carretero et al. 20151. We 
modify and extend this algorithm in two ways: (i) we corre¬ 
late galaxy colours with the parent halo concentration and 
(ii) we calculate a stellar mass for each galaxy using a colour- 
dependent mass-to-light ratio. 


2.2 Correlating colour and concentration 

As indicated above, in the basic S09 algorithm the colour of 
any given galaxy is assigned in two steps: first the galaxy 
is labelled ‘red’ (a ‘red flag’ is set to 1) with probability 
p(red|gal, Mr) where ‘gal’ is either ‘cen’ or ‘sat’, and then 
its colour is drawn as a Gaussian random number with the 
appropriate mean and variance. We correlate galaxy colours 
with parent halo concentrations by modifying the first step to 
adjust the red flag according to the concentration. The trend 
we wish to introduce is that more concentrated, older halos 
should host older, or redder, galaxies. 

At fixed halo mass, halo concentration c is approximately 
Lognormally distributed, i.e. Inc is approximately Gaussian 
distributed with constant scatter uinc — 0.141n(10) and a 
mean value (Inc) = In c(m,z) that depends on halo mass and 


redshift (see, e.g. Bullock et al. 2001 Wechsler et al. 20021. 


In the range of interest the mean value is well described by 
( Ludlow et al.|[2014 l 


c(m, z) = a u{nn, z) 


( 1 ) 


where v(Tn,z) = 5c{z)/a{m) is the dimensionless “peak 
height” of the halc|^ and the coefficients a and /3 depend on 
the choice of halo mass definition. We will construct catalogs 
using the m 200 b definition for which the Zehavi et al. ( 2011| 
HOD was calibrated (m 2 oob is the mass contained in the ra¬ 
dius i? 200 b where the spherically averaged density is 200 times 
the background matter density p(z)) and the corresponding 
c(m, z) relation. For this choice we have a = 7.7 and j3 = 0.4 
(see, e.g.. Appendix A of Paranja^|2014 1. 

Let p(red) be the value of the red fraction in the absence 
of any correlation with concentration. (For clarity we drop 
the explicit dependence on galaxy luminosity and type.) It is 
useful to define the quantity 


g _ ln(c/c) 

^In c 


( 2 ) 


whose distribution p{s) is Gaussian with zero mean and unit 
variance, with c(m) and crinc defined above. We compute a 
conditional red fraction p(red|s) that depends on the parent 
halo concentration by setting 


p(red|s) = (1 - p)p(red) -I- p0(s - Sred), 


( 3 ) 


^ See |Skibba| ) |2009| l for a comparison of central and satellite colours 
using mocks based on the S09 algorithm and in an earlier version 
of the Y07 catalog. 

^ Sc{z) and cr{m) are, respectively, the critical density for spherical 
collapse at redshift z and the r.m.s. of initial fluctuations smoothed 
on mass scale m, each extrapolated to 2 = 0 using linear theory. 


where the Heaviside function 0(a:) is unity for a; > 0 and zero 
otherwise. This says that galaxies in halos with concentrations 
s > Sred have an enhanced probability to be red as compared 
to the S09 model, while this probability is lowered by a factor 
(1 — p) compared to p(red) for galaxies in low-concentration 
halos. The dividing line Sred between high and low concentra¬ 
tions is defined in such a way that the average red fraction 
across all halos satisfies 

/ OO 

dsp(s)p(red|s) = p(red), (4) 

-OO 

which, using equation |3f, implies 

p(red) = p(s > Sred) = —° • (5) 

Equation ([^ gives a step-like dependence of the red frac¬ 
tion on concentration; in principle one could also imagine 
schemes where the red fraction increases with concentration 
in a continuous manner, but in this case it becomes more 
complicated to simultaneously ensure (p(red| s)) = p(red) and 
0 < p(red|s) < 1, due to the Gaussian integral in equation (|^. 

The parameter p lies between zero and unity and controls 
the strength of the correlation between the red fraction and 
concentration; we discuss its physical meaning below. Setting 
p = 0 gives us the uncorrelated case when the red fraction 
does not depend on parent halo concentration. Setting p = 1 
on the other hand corresponds to ‘complete correlation’ where 
all galaxies in high(low)-concentration halos are labelled red 
(blue). The intermediate case 0 < p < 1 clearly interpolates 
between these two extremes. 

Having set up our model, we proceed by assigning colours 
to the galaxies. The red flag of the galaxy is chosen by drawing 
a uniform random number u £ [0,1) and setting the red flag 
to unity if u < p(red|s) and zero otherwise. Next, we draw 
a Gaussian random number g — r from the appropriate red 
or blue distribution. In the basic S09 algorithm, the latter 
step also does not depend on parent halo concentration. In 
principle, we could once again use the values of s to introduce 
a further correlation between the actual g — r values and halo 
concentration, but in practice it turns out to be difficult to do 
this while preserving the global colour distributions at fixed 
luminosity (which we must do since these are constrained by 
data). We therefore stick to the S09 procedure at this stage 
and draw from the appropriate ‘red’/‘blue’ Gaussian without 
further correlation with concentration. 

Ideally, we would correlate galaxy colours with the con¬ 
centrations actually measured in the simulation. However, 
these concentrations are only approximately Lognormal, and 
deviations from the Lognormal shape render the second equal¬ 
ity in equation f5f above invalid. Gontinuing to use s^ed = 
\/2erfc~^(2p(red)) then leads to the overall central and satel¬ 
lite red fractions not being preserved, with (p(red|s)) 
p(red). We could correct for this by calibrating the exact shape 
of p(s) and numerically inverting the relation p(s > Sred) = 
p(red) to obtain Sred- Note, however, that what we are really 
interested in is not the actual shape of the distribution of con¬ 
centrations, but only their ranking in halos of fixed mass. It 
is therefore easier to simply “Gaussianize” the measured log- 
concentrations, proceeding as follows. 

(i) We first randomly draw a Gaussian variate s for each 
halo in the catalog and derive a Lognormal concentration c 
using equations 0 and 0. 
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(ii) Next, we bin the halos in 16 equi-log-spaced mass bin^ 
and consider the lists of measured and Lognormal concentra¬ 
tions of halos in each bin. We rank order and reassign the 
Lognormal concentrations in a given bin according to the mea¬ 
sured concentrations in that bin. In detail, if there are N halos 
in a bin, then the halo with the largest measured concen¬ 
tration is assigned the largest Lognormal concentration in 
this bin, with j — 1 , 2 ,..., N . 

(iii) These reassigned Lognormal values of c are then used 
to distribute satellites in their respective parent halos, and 
the corresponding Gaussian values of s are used to induce 
conformity in the colours as described above, with Sred given 
by equation 

This ensures that galaxy colours in this model are assigned 
consistently and inherit the large scale environment depen¬ 
dence of assembly bias. 

We will show that this procedure, which defines our de¬ 
fault model, leads to differences in the satellite red fraction 
in groups with blue and red centrals (1-halo conformity) as 
well as specific spatial trends of galaxy red fractions around 
blue and red isolated objects well outside i? 200 b (2-halo con¬ 
formity). As mentioned in the Introduction, it is also useful 
to ask how these trends would change if conformity arises not 
because of assembly bias, but due to some other halo property 
that does not show environmental trends. To this end, we will 
also consider a model in which galaxy colours are chosen and 
satellites are distributed exactly as in the default procedure 
above, except that we skip the rank ordering step (ii). This 
is equivalent to randomizing the halo concentrations at fixed 
halo mass before correlating them with the galaxy colours. In 
this model - which we will denote “no-^/i” - we expect our 
mocks to exhibit 1-halo conformity without a corresponding 
2 -halo signal. 

The parameter p has an interesting interpretation. It is 
easy to show that p satisfies the relations 

_ p(red|s > Sred) - p(red) 

^ i — p(red) 

^ p(red|old) -p(red) , , 

1 — p(red) 


where the second approximation holds if we assume that high 
concentration halos are old. Although simply a restatement of 
the definition of p in equation , equation ([^ allows us to 


interpret p as a “quenching efficiency” ( 

van den Bosch et al. 

2008 

|Peng et al.|20101120121 |Kovac et a 

.|2014| IPhillips et al. 

2014 

Knobel et al. 20141. Written like 

this, in our default 


model p is the fraction of blue (or star forming) galaxies that 
became red ^r quenched) because their respective parent ha¬ 
los grew oldj In the no-2h model, however, the connection 
between halo ages and galaxy colours is lost (since the former 
exhibit halo assembly bias at large scales while the latter will 
not), and the second relation in equation will not hold. 
In each model of conformity, however, p can be thought of as 
a “group quenching efficiency”, driven by the age of the par¬ 
ent halo of the group in the default model and by some local 


^ The number of bins was chosen after testing for convergence. 

® This can only be approximately correct, since, e.g., the present 
day population of blue satellites is not in general representative 
of the progenitors of all present day satellites. Our results do not 
depend on this interpretation, however, and one could simply treat 
p as a free parameter in the model. 


but otherwise unspecified property of the group in the no-2h 
model. 

Our scheme of generating concentration-dependent red 
flags smoothly interpolates between the uncorrelated case 
(p = 0) and the nearly ‘completely correlated’ case (p = 1). 
In practice, we perform this operation separately for satel¬ 
lites and centrals, which ensures that the overall central and 
satellite red fractions are left unaltered by construction. The 
model is flexible enough that, if needed, the group quenching 
efficiency p can be set separately for centrals and satellites, 
with the interpretation that different quenching mechanisms 
might play a role for centrals and satellites in the same halo. 
E.g., this might be the case if srtfchalo age is more relevant 


for satellite colour than the age of the parent halo (Hearin & 
|Watson||2013p . In this work, however, we do not explore this 
possibility and only use a one-parameter setup with p taken to 
be the same for centrals and satellites, independent of galaxy 
luminosity and halo mass. Galactic conformity arises from the 
fact that the colours of centrals and satellites in a group are 
affected by their common group concentration. The fact that 
the concentrations have a scatter at fixed halo mass means 
that the conformity will persist even when binning objects by 
the parent halo mass. 


2.3 Calculating stellar masses 

As a separate step, having assigned galaxy colours, we com¬ 
pute stellar masses using a colour dependent mass-to-light ra¬ 
tio that we have also fit to SDSS galaxies. In particular, we 
use 


{M/L)r = 4.66(5 ~ C ~ 1-36(5 “ ~ L108, 


( 7 ) 


to compute the °'^r-band Petrosian mass-to-light ratio 
{M/L)r in units of Mq/Lq, and add a Gaussian 1-sigma scat¬ 
ter of 0.2 (we describe the procedure and data set used to ob¬ 
tain this fit in section]^. Taking the base-10 logarithm and 
adding {Mr — 4.76)/(—2.5) gives us the log stellar mass in 
units of h~‘^MQ. This does not guarantee that the central of a 
group, which is the brightest by construction, is also the most 
massive. In practice this affects about 8-10% of the groups in 
the total sample, with the actual number depending on the 
the value of p. This may be related to the findings of |Skibba| 
et al. (20111 regarding a non-zero fraction of groups in which 


the central is not the brightest; we will explore this further in 
future work. 


3 SIMULATIONS 


Our mocks are built on a suite of A-body simulations of cold 
dark matter (GDM), each of which resolves a periodic cubic 
box of comoving volume (200/i~^Mpc)® with 512^ particles 
using our fiducial cosmology. This gives us a particle mass of 
nipart = 4.1 X 10^h~^M q. The simulations were run using the 
tree-PM cod^ Gadget-2 (Springel 20051 with a force reso¬ 
lution £ = 12.5/i“^kpc comoving (~ 1/30 of the mean particle 
separation) and a 1024® PM grid. Initial conditions were gen¬ 
erated a t z = 99 employing 2”‘*-order Lagra^ian Pert urbation 
Theory ( Scoccimarro|| 19981, using the cod^ Music (Hahn & 


^ http://www.mpa-garching.mpg.de/gadget/ 

^ http://www.phys.ethz.cIi/~hahn/MUyUJ/ 
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Abel||20rT I with a transfer function calcnlated nsing the pre¬ 


scription of Eisenstein & Hu (19981. We have run 10 realisa¬ 
tions of this simulation by changing the random number seed 
used for generating the initial conditions and have used the re¬ 
sulting 2 = 0 snapshots for our mocks. The simulations were 
run on the Brutus clustei[3at ETH Ziirich. 

To identify halos, we have used the cod^ Rockstar 
( Behroozi et al.|2013b (, which assigns particles to halos based 
on an adaptive hierarchical Friends-of-Friends algorithm in 
6-dimensional phase space. Rockstar has been shown to 
be robust for a variety of diagnostics such as density pro¬ 
files, velocity dispersions, merger histories and the halo mass 
function. As mentioned earlier, we use the m 2 oob value re¬ 
ported by Rockstar as the mass of the parent halo, and 
Rzoob as its radius. We use the value of the halo scale ra¬ 
dius rg reported by Rockstar to compute halo concentration 
C 2 oob = R 2 oob/Ts. The smallest halo mass we resolve sets the 
faintest luminosity we can reliably sample. We discard objects 
having m 2 oob < 20mpart, which allows us to set the lumi¬ 
nosity threshold to Mr,max = —19.0. Having 10 independent 
simulations means that we can reliably characterise the sample 
variance on our observables for each mock configuration. 


4 DATA 

In this section we describe the data sets we use for comparison 
with the mocks, and give details of various fits used in dehning 
the mocks. 


4.1 The YOT catalog 

We base some of the fits required for generating our mocks (the 
colour distribution and satellite red fraction at fixed luminos¬ 
ity, and the mass-to-light ratio at fixed colour) on the galaxies 
contained in the Y07 group catalo j^ These fits do not depend 
on the level of conformity, and we therefore additionally use 
the Y07 catalog as a baseline for galactic conformity as well. 
The Y07 catalog is built using the halo-based group finder de¬ 
scribed in [Yang et al.| ( |2005| ) to identify groups in the New 
York Universi ty Val ue Added Galaxy Cata log (NYU-V AGC; 
Blanton et al.|[2005 (, based on the SD S^^ ( ^rk et al.||2000 l 
data release 7 (DR7; [Abazajian et al.|[2009[ ). Throughout, we 
will only consider galaxies with spectroscopic redshifts (‘sam¬ 
ple IF of Y07) restricted to the range 0.01 < z < 0.07, and 
also restrict the Petrosian absolute magnitudes to the range 
-23.7 <Mr < -16. 

In part of the analysis below, we will use a luminosity- 
complete subsample - denoted “Y07-Lum” - containing galax¬ 
ies with Petrosian absolute magnitudes Mr < Mr,max = —19.0 
and Model g — r colours from the NYU-VAGC as provided by 
Y07. Additionally, when comparing properties as a function of 
stellar mass, we will use another subsample - denoted “Y07- 
Mass” - containing values of Petrosian stellar mass m* and 
Model g — r colours obtained by ru nning the kcorrect_v4.2 
cod of [Blanton fc Roweis (20071 on the corresponding Pet¬ 
rosian and Model properties, respectively, of the Y07 galax- 


® http://www.cluster.ethz.ch/index_EN 
® http://code.google.com/p/rockstar/ 

http://gax.shao.ac.cn/data/Group.html 
http://www.sdss.org 

http://howdy.physics.nyu.edu/index.php/Kcorrect 






Figure 1. Comparison of observed Model g — r distributions in 
bins of Petrosian Mr in the Y07-Lum catalog (solid green his¬ 
tograms) with the corresponding double-Gaussian fits reported in 
equations §. We show the ‘red’ and ‘blue’ components in the fits 
as the respectively red and blue dashed lines, and the total of these 
as the dashed green lines. For comparison we show the S09 fits as 
dotted lines; the discrepancy between S09 and the Y07 catalog is 
most likely due to differences between DR7 and DR4 (which was 
used by S09) and our choice of using Model colours. 


ies obtained from the MPA-JHU cataloJ^ 


(Kauffmann et al. 


20031 by matching galaxy IDs across the MPA-JHU and NYU 


catalogs. We use magnitudes that are K-corrected to their 
values at t: = 0.1. We have checked that there is reasonable 
agreement between comparable galaxy properties in the NYU- 
VAGC and the MPA-JHU catalog. When showing mass func¬ 
tions and fractions, we apply inverse Vjnax-weighting to the 
galaxies in the Y07-Mass catalog. 

Our choice of Petrosian quantities for defining both the 
stellar masses as well as the mass-to-light fit (see below) in the 
Y07-Mass catalog is motivated by the fact that the HOD we 
use has been calibrated on Petrosian absolute magnitudes by 


Zehavi et al. (20111. We have checked that there is reasonable 


agreement between Model and Petrosian stellar masses in the 
Y07-Mass catalog and, moreover, that the trends in the red 
fractions and associated observables which we discuss below 


http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Figure 2. Conformity-independent variables as a function of luminosity (left panel) and stellar mass (right panel). Filled symbols joined by 
solid lines show the results of averaging over 10 realisations of our mocks, while open symbols joined by dashed lines show measurements 
in the Y07 catalog. The error bars on the mock results show the r.m.s. fluctuations (standard deviation) around the mean value, while the 
errors on the Y07 data are estimated from 150 bootstrap resamplings. The cyan diamonds show the satellite fraction /sat- The dark red 
crosses, green circles and yellow triangles respectively show the red fractions for all galaxies (/red)? only centrals (/red|cen) nnd only satellites 
(/red|sat)- The magenta squares show the satellite quenching efficiency ^sat = (/red|sat ~ /red|cen)/(l ~ /redlcen)- The ‘red’/‘blue’ classification 
for both data and mocks uses equation ( |10[ | for the left panel and equation for the right panel. 


differ by only ~ 10% when using Model or Petrosian quan¬ 
tities. Consequently, we do not expect any of our qualitative 
conclusions to be affected by our choice. 


4.2 Mass-to-light ratio fit 

We fit a quadratic relation to the ®'^r-band Petrosian mass- 
to-light ratio {M/L)r at fixed value of Model g — r colour in 
the Y07-Mass catalog. In particular, we compute the median 
value of {M/L)r in bins of g — r, using the median g — r as 
the bin center. The resulting least-squares best fit is given in 
equation 0 and is close to but slightly lower than the one in 
equation (1) from Wang & White 120121. The measured r.m.s. 
scatter of {M/L)r in the bins of g — r has an average value of 
0.2 which we also use as described below equation 0 . 


4.3 Color-luminosity fits 

We fit double-Gaussian shapes to the distributions of Model 
g — r colours in bins of Petrosian absolute magnitude Mr in the 
Y07-Lum catalog. The resulting fits can be summarised using 
5 quantities: the means and variances of the ‘red’ and ‘blue’ 
distributions and the probability p(red|Mr) that the colour is 
drawn from the ‘red’ distribution. The full g — r distribution 
can then be written as 


for the blue distribution. We find the following best fit values: 

piredlMr) = 0.423 - 0.175 (Mr + 19.5) 

{g-r\Mr = 0.9050 - 0.0257 {Mr + 19.5) 
{g-r\Mr = 0.575 - 0.126 {Mr + 19.5) 
ar^d{Mr) = 0.0519 -f 0.0085 {Mr + 19.5) 
crbiue(M,.) = 0.150 + 0.015 {Mr + 19.5) (9) 

Figure[^compares these analytic functions (dashed lines) with 
the measured g — r distribution in the Y07-Lum catalog (solid 
histograms). For comparison, the dotted lines show the fits re¬ 
ported by S09 which were based on a similar catalog using DR4 
of SDSS; these are rather different from the Y07-Lum data, 
which could be due to differences between the DR7 colours 
and those in DR4. (The discrepancy reduces but does not dis¬ 
appear if we use Petrosian colours.) Overall, we see a good 
agreement between the data and the fits, except for a “green 
valley” which is not captured by the sum of two Gaussians. 
In principle, one can obtain a better fit by introducing a third 
‘green’ component ( [Krause et al.|20T^ [Carretero et al.j2015[ ), 
but we choose not to do this here. 

Although the algorithm classifies objects as red or blue 
depending on which Gaussian distribution their colours are 
drawn from, for ease of comparison with observational results, 
we will use sharp thresholds in g — r. When showing results as 
a function of luminosity we classify objects as red when their 
g — r colour exceeds 

(<? - r)cut = 0.8 - 0.03 (M,-b 20) , (10) 


p{g - r\Mr) = p{red\Mr)pTed{g - r\Mr) 

+ {I - p{red\Mr)) phiue{g - r\Mr), (8) 

where, e.g., Pred{g — r\Mr) is Gaussian with mean 
{g — r\Mr and standard deviation ared{Mr), and similarly 


and as blue otherwise ( jZehavi et al.|2005| ). When showing re¬ 
sults as a function of stellar mass we instead use the threshold 

{g - r)cut = 0.76 -b 0.10 [logio(m*) - 10] , (11) 

which is somewhat shallower than the threshold quoted by 
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Figure 3. All-galaxy luminosity (left panel) and stellar mass functions (right panel). Red squares show the average over 10 independent mock 
catalogs, with error bars showing the r.m.s. scatter. Cyan circles show measurements in the Y07-Lum (left) and Y07-Mass catalogs (right), 
with errors estimated from 150 bootstrap resamplings. For compar ison, the short-dashed black curve in the left panel shows the Schechter 
fit to the '^■^r-band luminosity function in SDSS data reported by [Blanton et ah (2003 i, while the corresponding curve in the right panel 

The vertical dotted line in the right panel shows 


shows the stellar mass function fit for the Y07 catalog as reported by [Peng et al. 
the approximate mass-completeness limit of our mocks at log]^Q(m*) = 9.9. 


(2012 


van den Bosch et al. (20081 who use a slope of 0.15 instead 


of 0.1. We have found that equation (111 gives a slightly bet¬ 
ter separation of the bi-modal colour-mass distribution in the 
Y07-Mass catalog than the van den Bosch et al. (2008) rela¬ 
tion does, although our qualitative conclusions do not depend 
on the choice of threshold. We apply the same threshold to 
colours in our mocks as well as in the Y07 data. Using sharp 
thresholds means that the measured red fractions /red will, in 
general, differ from the probabilities p(red) discussed earlier. 


5 RESULTS 

We now present the results of our mock algorithm and compare 
them with corresponding measurements in the Y07 catalog. 
We start by discussing observables that do not depend on the 
presence or absence of conformity. 


5.1 Conformity-independent observables 

The following observables are independent of the level of con¬ 
formity (i.e., the value of p) because their definition does not 
involve a simultaneous determination of central and satellite 
colour: the fraction of galaxies that are satellites (/sat), the 
average red fractions of all galaxies, centrals and satellites (re¬ 
spectively, /red, /red|cen and /redisat), and the Satellite quench¬ 
ing efficiency Seat = (/red|sat “ /red|cen)/(l “ /red|cen)- AmOUg 
these functions, the luminosity dependence of /sat is deter¬ 
mined solely by the HOD (which itself is fit using the lumi¬ 
nosity dependence of clustering), that of /red is fixed by the 
double-Gaussian fitj^ and the remaining quantities addition- 


Note that /red depends on the actual shape of the double- 
Gaussian, not only on p(red|Mr.). 


ally depend on the choice of satellite red fraction which we 
discuss next. 

As mentioned previously, we restrict ourselves to a simple 
prescription in which p(red|sat, Mr) (the probability that g — r 
for a satellite of luminosity Mr is drawn from the ‘red’ Gaus¬ 
sian) as well as the corresponding quantity p(red|cen, Mr) for 
centrals are independent of halo mass, in which case only one 
of these can be chosen independently. By comparing to the 
measured luminosity dependence of /red | sat, we find reason¬ 
able agreement using 


p(red|sat, Mr) = 1.0 — 0.33 


1 + tanh 


/Mr+ 20.25 

V 2T 


, ( 12 ) 


although this is by no means the only function that can do 
the job. This in turn fixes the luminosity dependece of /red|cen 
and Esat- The stellar mass dependence of all the quantities 
mentioned above is then also completely fixed. 

The results of the comparison are shown in Figure]^ as 
a function of luminosity in the left panel and stellar mass in 
the right panel. The open symbols joined by dashed lines show 
the measurements in the Y07-Lum (left panel) and Y07-Mass 
catalogs (right panel), while the filled symbols joined by solid 
lines show the mean over 10 realisations of our mocks. The 
error bars on the mock data show the r.m.s. scatter (sam¬ 
ple variance) of the respective quantities over all realisations, 
while the error bars on the Y07 data are estimated from 150 
bootstrap resamplings. We see that, apart from systematic dif¬ 
ferences of < 10%, there is good agreement between the mocks 
and the data. In particular, the mocks correctly reproduce the 
observed near-independence of Esat on stellar mass and lumi¬ 
nosity, as well as matching its amplitude. The systematic dif¬ 
ferences could be due to one or more of the following: residual 
systematics in the HOD parameters and the mass-to-light fit, 
the fact that we do not exactly reproduce the observed colour 
distribution, and our choice of halo mass-independent central 
and satellite red fractions. In particular, the difference in f^ed 
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Figure 4. Number counts of centrals as a function of g — r and 
luminosity (top panel) or stellar mass (bottom panel) in one of our 
mocks. We clearly see how a luminosity threshold propagates into a 
mass incompleteness for red objects. A visual inspection indicates 
that our mocks should be mass complete for log]^Q(m*) > 9.9 (ver¬ 
tical dotted line in the bottom panel). The thick yellow lines show 
equation in the top panel and equation ED in the bottom 

panel. For comparison, the dashed yellow line in the bottom panel 
shows the threshold used by [van den Bosch et al.| ||2008[l. 


is very likely due to the mismatch in the colour distributions 
around the “green valley” (see Section 4.3 and Figure[^. Since 
we match the satellite red fraction /red | sat quite well by con¬ 
struction, the difference in /red then also causes a difference 

in /red|cen- 

Finally, an important set of observables that are inde¬ 
pendent of conformity are the luminosity and stellar mass 
functions of the mocks. Figure compares these with cor¬ 
responding measurements in, respectively, the Y07-Lum and 
Y07-Mass catalogs. We also show the fits to SDSS data re¬ 


ported by Blanton et al. (20031 for the luminosity function 
and Peng et al. (20121 for the mass function. Overall there 
is good agreement between the measurements in our mocks, 
in the data and the corresponding fits from previous work; 
this serves as a sanity check on our algoritm. The fact that 
our measurements in the Y07-Mass catalog do not quite agree 
with the fits reported by Peng et al. (20121 highlights the fact 
that the shape and amplitude of the mass function is quite 
sensitive to the exact definition of stellar mass. 

From a visual inspection, the mocks appear to be mass 
complete for log]^Q(m*) > 9.9, which is the threshold we will 
use below. For our choice of h = 0.7 in the Y-body simulations, 
this corresponds to logjQ(m*/M 0 ) > 10.2. The incompleteness 
at lower masses is because our mock is luminosity complete 
with Mr < —19.0; the colour-dependence of the mass-to-light 
ratio ([^ then means that we are missing most of the faint 
red galaxies (and many of the faint blue ones) that would 
populate the low mass end. This becomes clearer in Figure]^ 
which shows number counts of centrals as a function of gr — r 
and luminosity (top panel) or stellar mass (bottom panel) in 
one of our mocks. The deficit of low-mass red objects is clear in 
the bottom panel. Similar results hold for the satellites as well. 
Our choice of mass-completeness threshold is further justified 
by the fact that only ~ 1 % of the galaxies with logjQ(m*) > 
9.9 belong to the faintest bin —19.2 < Mr < —19.0, implying 
that fainter galaxies would contribute negligibly above this 
mass threshold. 


5.2 Effects of correlating galaxy colour and halo 
concentration 

We now investigate the effects of a non-zero value of p in ob¬ 
servables that do respond to a correlation between central and 
satellite colours. To start with we simply explore the effects of 
changing the value of p, and later motivate a specific value by 
comparing to the Y07 catalog. 

The left panel of Figure shows /red|sat as a function of 
satellite stellar mass, split by whether the corresponding cen¬ 
tral is red or blue, for two choices of p = 0.01,0.99. For each 
of these choices, the points show the mean over 10 indepen¬ 
dent mocks. Since we are interested in understanding these 
mean trends in this plot, the error bars in both panels show 
the standard error on the mean (rather than the r.m.s. scat¬ 
ter). The plot clearly shows that a positive correlation goes in 
the direction of explaining conformity, which has been noted 
earlier by others as well. The open symbols (red triangles and 
blue inverted triangles) joined by dotted lines show the result 
of the algorithm setting p = 0 . 01 , i.e. essentially no corre¬ 
lation, in keeping with what is usually assumed by standard 
HOD based algorithms. In this case there is no distinction 
between the two red fractions at the bright end, but fainter 
satellites with blue centrals are preferentially slightly redder 
than those with red centrals (a weak ‘negative conformity’). 
We have checked that the latter effect is due to averaging over 
luminosities and completely disappears if we plot results at 
fixed luminosity insteacf^ (not shown). 

The filled symbols with solid/dashed lines show the result 


This traces back to our choice of satellite red fraction which is 
independent of halo mass. Had this not been the case, we would have 
seen a similar effect in the luminosity-based plot as well, arising in 
this case from an averaging over halo mass. The effect in the stellar 
mass plot would now be even more pronounced. 
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Figure 5. (Left panel): Satellite red fraction as a function of satellite stellar mass, split according to whether the associated centrals are red 
(red points and lines) or blue (blue points and lines). The open red triangles and blue inverted triangles with dotted lines show the result 
of the algorithm setting p = 0.01, i.e. essentially no correlation between galaxy colour and host halo concentration, as is routinely assumed 
by standard HOD based algorithms. The filled symbols with solid/dashed lines show the result when p = 0.99, i.e. strong correlation, in the 
presence (default model; red triangles and blue inverted triangles with solid lines) or absence (no-2h model; red crosses and blue stars with 
dashed lines) of 2-halo conformity as described in the text. The points show the mean over 10 independent mocks and the error bars the 
standard error on the mean. (Right panel): The red fraction of all galaxies surrounding red or blue “isolated primaries” as defined in the 
text, as a function of spherical distance, colour-coded and formatted as in the left panel, for the same set of 10 mocks used for the respective 
measurements in the left panel. All objects are required to have logiQ(m*) > 9.9 and we do not restrict the masses of their parent halos. The 
horizontal cyan line and band respectively show the average all-galaxy red fraction above this stellar mass threshold and its r.m.s. scatter 
over 10 independent mocks (these values are independent of p). Note that the scale on the vertical axis differs from that in the left panel. 
For both centrals and satellites, the distinction between red and blue was made using equation 



Figure 6. Same as Figure]^ except that the filled red crosses and blue stars joined by dashed lines now show the result for mocks with 
the default model, but where we set p = 0.5 (“medium” correlation between colour and concentration). The open triangles with dotted lines 
and filled triangles with solid lines are the same as in Figure]^ This plot demonstrates the tunability of our model, which is one of our key 
results. 
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when p — 0.99, i.e. strong correlation. In this case satellites 
with red centrals are clearly significantly redder than satellites 
with blue centrals. The filled red triangles and blue inverted 
triangles joined by solid lines show the result for mocks that 
used our default model of conformity, rank ordering the Log¬ 
normal halo concentrations according to the measured con¬ 
centrations from the simulation in 16 equi-log-spaced mass 
bins (see Section 2.2 1 . The red crosses and blue stars with 
dashed lines show the result for mocks that used the no-2h 
model, with randomly assigned Lognormal concentrations. As 
expected, the red fractions for these two cases are nearly iden¬ 
tical. 

The right panel of the Figure is formatted identically to 
the left panel and uses the same set of mocks. In this case, we 
show the red fraction of all galaxies with logjQ(m*) > 9.9 as 
a function of their distance r from red or blue “isolated pri¬ 
maries”; we define the latter as galaxies of mass m* that do 
not have any galaxy more massive than m*/2 within a spher¬ 
ical radius of SOOkpcf^ We find that the galaxies picked by 
this definition are predominantly centrals, with ~ 10% of iso¬ 
lated primaries being satellites (see also Hearin et al. 20151. 
We consider isolated primaries in order to be close to the se¬ 
lection criteria used in recent observational studies (see, e.g.. 


Kauffmann et al. 2013We investigate the impact of changing 
the selection criterion in Appendix |A2[ 


The horizontal cyan line and associated band respectively 
show the mean and r.m.s. scatter over 10 realisations of the 
average all-galaxy red fraction which is independent of con¬ 
formity strength. The mocks with p = 0.01 show nearly iden¬ 
tical trends at r > 2Mpc regardless of the colour of the iso¬ 
lated primary, although the red fractions remain above the 
global value even as far as lOMpc from the center. There is 
a weak, conformity-like signal at small separations, which we 
have checked is entirely due to averaging over halo mass and 
disappears when using galaxies in fixed bins of halo mass. To¬ 
gether with the apparent ‘negative conformity’ seen for this set 
of mocks in the left panel, this reiterates the need to be cau¬ 
tious when interpreting trends in analyses that average over 
luminosity and halo mass. 

Unlike the left panel, the two sets of mocks with p — 0.99 
now show distinct trends. In the no-2h mocks, the red fractions 
around blue and red isolated primaries are different until a dis¬ 
tance of ~ 4Mpc from the center (the one around red isolated 
primaries being higher), beyond which they are identical and 
also coincide with the “no conformity” red fractions associated 
with p = 0.01. The default mocks, on the other hand, show 
a similar but substantially larger difference between the red 
fractions out to ~ 6Mpc, beyond which they are also nearly 
identical but substantially larger than the red fractions in the 
other data sets (see sectionj^and Appendix A2 for a discussion 
of these trends). 

Figure shows the same quantities as Figure except 
that the filled red crosses and blue stars joined by dashed 
lines now show the results for default mocks where we set 
p = 0.5 (medium correlation). This plot shows how the dif- 


This definition is somewhat different from the observationally- 
motivated one employed by [Kauffmann et ah] | |2013| l , who used a 
cylinder in redshift space and projected distance (see also |Hearin| 
jet al.|2015| . In this paper, however, we do not compare our 2-halo 
results with observations, and the simpler 3-dimensional definition 
above suffices to understand various trends in the signal. The signal 
itself is expected to diminish in strength due to projection effects. 



Figure 7. Comparison of “conformity red fractions” /red|sat re 
/redlsat.bc of satellites having red and blue centrals, respectively, 
with the red fraction of centrals /redlcem for the Y07-Mass catalog 
(open symbols joined by dashed lines) and in mocks that used p = 
0.65 (filled symbols with solid lines). The red triangles, blue inverted 
triangles and green circles respectively show /re^sat,rc, /red|sat,bc 
and /red|cen Same as in Figure pi. The error bars 

on the mock results show the r.m.s. fluctuations around the mean 
value over 10 realisations, while the error bars on the Y07 data are 
estimated from 150 bootstrap resamplings. Notice the remarkable 
similarity between /red|aat,bc s-'id /red|ceni especially at small masses 
in the Y07-Mass catalog. 


ference between the red fractions around bine and red cen¬ 
trals/isolated primaries changes with group quenching effi¬ 
ciency p and demonstrates the tunability of our model, which 
is one of our key results. 


5.3 Fixing the level of conformity 

We notice in the left panel of Figure that the red fraction 
/red|aat,bc of Satellites with blue centrals is particularly sensi¬ 
tive to the value of p. This happens, at least in part, due to 
the inherent asymmetry of satellite colours, most of which are 
red. Consequently, we can use measurements of 1-halo confor¬ 
mity in the Y07 catalog to hx the value of p, which can then 
be used to predict the 2-halo conformity signal. 

By trial and error we have found that setting p within 
~ 10% of p = 0.65 leads to a behaviour of /red|sat,bc in the 
mocks that closely resembles the one in the data over the 
mass range allowed by the mocks. Figureshows red fractions 
averaged over 10 mocks that used p = 0.65, and we see that 
both /red|sat,bc and /redisat.rc agree very well with the Y07 
measurements. The results of the previous section show that 
this agreement is independent of whether we use the default or 
the no-2h model. Notice that the red fraction of satellites with 
blue centrals /red|sat,bc in both data and mocks is remarkably 
similar to the corresponding /red|cen. This is particularly true 
for the data at low masses where our current mocks do not 
reach. 

Figure shows the large scale conformity signal predicted 
in the mocks having p = 0.65, using both the default and the 
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Figure 8. Large scale conformity predicted in three configurations 
of mocks that used p = 0.65 and whose respective 1-halo conformity 
signals are each consistent with measurements in the Y07 catalog. 
Filled symbols (red triangles and blue inverted triangles) with solid 
lines used our default mocks with halo mass-independent satellite 
and central red fractions at fixed luminosity, in which 2-halo con¬ 
formity is switched on. Filled symbols (red crosses and blue stars) 
with dashed lines used the no-2h mocks in which 2-halo confor¬ 
mity is switched off. Finally, open symbols (red triangles and blue 
inverted triangles) with dotted lines used mocks in which 2-halo 
conformity was switched on and the red fraction of centrals at fixed 
luminosity had an additional dependence on halo mass as described 
in Appendix |Al| The error bars on the mock results show the r.m.s. 
fluctuations around the mean value over 10 realisations for each con¬ 
figuration. The inset zooms in on the behaviour of the signal in the 
default mocks at very large separations. Unlike Figures and 
we restrict the isolated primaries to lie in the stellar mass range 
10.0 < log]^Q(m*) < 10.5, while their neighbours are required to 
have log]^Q(m*) > 9.9. The horizontal cyan line and band respec¬ 
tively show the average all-galaxy red fraction above the latter mass 
threshold and its r.m.s. scatter over 10 independent mocks (these 
values are independent of p). 


no-2h models. To assess the effect of our assumption of mass- 
independent satellite and central red fractions, we also show 
the corresponding signal in mocks having 2-halo conformity, 
but in which the red fraction of centrals is mass-dependent as 
described in Appendix |A1| These were constructed to match 
the conformity-independent observables discussed earlier, as 
well as the 1-halo conformity signal, for which it suffices to 
use p — 0.65 again (see Figure All. In each case, we restrict 
the stellar mass of the isolated primaries to lie in the range 
10.0 < logjQ(m*) < 10.5. We see that the overall signal in all 
three mocks is very similar at scales < 4Mpc. At larger scales, 
the mocks with a halo mass dependent central red fraction 
show a slightly higher signal than in our default mocks. The 


large scale signal in the no-2h mocks, on the other hand, is 
noticeably smaller than in the default. The inset zooms in 
on the behaviour of the signal in the default mocks at very 
large separations. The error bars depict the r.m.s. fluctuations 
around the mean over 10 realisations for each configuration. 

Figure shows a break-down of these large scale trends 
as a function of halo mass. Each panel shows the results for 
our default and no-2h mocks (formatted identically to Figure]^ 
and with the same restrictions on stellar masses of the isolated 
primaries and their neighbours), with the top left panel show¬ 
ing the same configuration as Figure]^ and the top right, bot¬ 
tom left and bottom right panels showing bins of increasingly 
larger halo mass. Each panel shows the average all-galaxy red 
fraction in the respective halo mass bin, and its r.m.s. scat¬ 
ter, as the horizontal cyan line and band, respectively. We 
clearly see that the large difference between red fractions at 
separations < 4Mpc arises primarily from large halo masses, 
while the distinction between the default and no-2h signal at 
> 8Mpc only emerges at the smallest halo masses. We discuss 
these results in section [6] below. 


6 DISCUSSION 


The 1-halo conformity trends seen in the left panels of Fig¬ 
ures [^and[§ and in the mock results in Figurej^are a straight¬ 
forward consequence of our choice of implementation and the 
definition of the group quenching efficiency p as discussed in 
section As p increases and approaches unity in our model, 
red (blue) satellites preferentially live in halos hosting red 
(blue) centrals. An interesting aspect of Figure is the re¬ 
markable similarity of the red fraction of satellites with blue 
centrals and the red fraction of centrals at low stellar mass. A 


similar result was obtained by Phillips et al. (20141, who found 


that the SFRs of satellites with star forming centrals are essen¬ 
tially the same as those of isolated galaxies with similar stellar 
mass. This would imply that the quenching efficiencies SsatiSFc 
of the satellites with star forming centrals are consistent with 


being zero. On the other hand, Knobel et al. (20141 reported 


values of SsatiSFc different from zero. This apparent discrep¬ 
ancy is possibly due to differences in the sample definitions in 


these two studies; in particular, unlike Phillips et al. (20141, 
[Knobel et al.j ( |2014[ ) limited their analysis to satellites residing 
in groups with at least 3 members with m* > IO^^Mq, thereby 
removing many low mass groups from the sample. As [Knobel] 
jet al.j p5T4| ) also report a strong positive dependence of £aat 
on halo mass, it is natural to expect that the quenching effi¬ 
ciency in a sample biased toward more massive haloes would 
be nonzero. It will be very interesting to extend our mocks to 
lower stellar masses using a higher resolution simulation and 
compare the resulting 1-halo conformity trends - particularly 
those at fixed group richness and halo mass - with the data, 
but this is beyond the scope of the present work. 

Turning to the large scale trends of red fractions surround¬ 
ing isolated primaries, we encounter an even richer picture. It 
has been argued in the recent literature that large scale differ¬ 
ences between red fractions or mean star formation rates sur¬ 
rounding blue and red isolated primaries must arise from envi¬ 


ronmental correlations across different halos (Kauffmann et al. 


20131 and are therefore a signal of 2-halo conformity (Hearin 
et al.|[^15p. Our results in Figure]^ and especially Figure [^ 


however, show that an alternative explanation is also possible. 
Figure [^ used three different sets of mocks, all of which are 
constructed to give a 1-halo conformity signal consistent with 
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Figure 9. Large scale conformity signal for p = 0.65, broken down by halo mass as indicated. Note that the mass restrictions apply to the 
parent halos of both the isolated primary as well as its neighbours. Each panel is formatted identically to Figure]^ and shows the results for 
the default and no-2h mocks. See text for a discussion. 


measurements in the Y07 catalog. Only two of these mocks 
have a genuine 2-halo conformity built in, however; the third 
contains randomly assigned halo concentrations which are not 
correlated with the halo environment. This third set of no- 
2h mocks serves as a useful toy example of the possibility that 
conformity might arise due to some property of halos that does 
not exhibit large scale environmental effects (see the discussion 
in the Introduction). The large scale conformity signal for all 
three mocks is very similar at scales < 4Mpc. Figure]^ shows 
that the difference of red fractions at these scales is dominated 
by halo masses logig(m) > 13.25 (bottom right panel) - for 
which it is consistent with being a 1-halo effeclp^ - and has a 
smaller and noisier contribution from small halo masses (top 
right panel, see also Appendix A2|. This means that, at these 
scales, we could easily mistake a purely 1-halo effect for gen¬ 
uine 2-halo conformit}!^ Similar conclusions can be drawn by 
comparing the p = 0.99 measurements in the right panel of 
Figure with those in the mocks with p = 0.01 or no con¬ 
formity. The p = 0.99 no-2h measurements are identical to 


For comparison, a cluster-sized halo of mass Mq has a 

radius i? 2 oob — l.TMpc for our cosmology, and the largest halos in 
our simulations reach masses of ~ 

Figure|^also shows that the large scale conformity signal appears 
to be robust against a mass dependence of the central red fraction, 
showing a slight increase on average as compared to the default 
mocks. 


the p = 0.01 measurements at large scales, which is expected 
since neither of these mocks has any 2-halo conformity. The 
measurements in the two p — 0.99 mocks on the other hand 
are very similar at scales < 4Mpc. Taken together, this shows 
that that the conformity signal at < 4Mpc in the p = 0.99 
mocks is almost entirely a 1-halo effect (see Appendix A2 for 
further discussion). 

From the observational point of view, halo mass- 
dependencies can in principle be controlled by performing tests 
at, e.g., fixed group richness ||Koester et al.|2007[|Skibba et ah 


2007 


Skibba 20091, stellar mass of the central (More et al. 


20111 or any other observable property that correlates well 
with halo mass (|Lu et al.|2015|). Notice, however, that all the 


measurements in Figures and use isolated primaries with 
stellar masses 10.0 < logjQ(m*) < 10.5 and still cannot distin¬ 
guish between the default and no-2h models at scales < 4Mpc, 
unless the halo mass is explicitly restricted to small values. We 
have also checked that the relative behaviour of our models at 
these scales remains unchanged when using primaries of dif¬ 
ferent stellar masj^ This is unfortunate, since it very likely 
means that one cannot conclusively argue that the confor¬ 
mity measured by|Kauffmann et al.|(|201^ at projected scales 


Figure [A3| in the Appendix shows that the signal, particularly 
at large scales, is essentially independent of the stellar mass of the 
primary; we attribute this to a large scatter in the distribution of 
halo masses at fixed stellar mass (Figure A4l. 
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< 4Mpc in similar mass bins is evidence of galaxy assembly 
bias. Projection effects are unlikely to alter this conclusion. 
It will be useful to repeat such tests at fixed group richness, 
which might be a better indicator of halo mass. 

For now, it is very interesting to note that there is, in fact, 
a signal which diminishes in the absence of genuine 2-halo con¬ 
formity and is also robust against the inclusion of halo mass 
dependence in the red fractions. This is the fact that, in the 
presence of 2-halo conformity, the red fractions at even larger 
scales (> 8Mpc) surrounding both red and blue isolated ob¬ 
jects remain above the global average value by a small but 
statistically significant amount, showing nearly identical val¬ 
ues and a trend that is relatively insensitive to the halo mass 
dependence of the central red fraction. When genuine 2-halo 
conformity is absent, the red fractions at these scales are dis- 
cernably closer to the global average; as noted earlier, they 
are also identical to the “no conformity” red fractions at these 
scales. These trends can be seen in Figures [Sl [8| and |A2[ with 
the elevation compared to the global value seen clearly in the 
zoom-in inset panel of Figure]^ Figurej^shows that this signal 
is absent at large halo masses and only emerges in the range 
11.25 < logio(m) < 12.25. 

This dependence on halo mass is exactly what one expects 
from a signal driven by halo assembly bias; in fact, one can 
also derive a more detailed understanding of the nature of 
the signal using the Halo Model. In Appendix |A2| e.g., we 
argue that the red fractions at very large separations from 
central galaxies are determined by ratios of cross-correlation 
functions of red/blue galaxies with red/blue centrals, in which 
the concentration-dependence of halo assembly bias couples 
with that of galaxy colour introduced in our default model 
and leads to an elevation that is qualitatively similar to what 
we see in Figure The fact that the signal is only of order a 
few percent is consistent with halo assembly bias being a weak 
effect. 


Our model assumes that central galaxies are always the 
brightest members of their respective groups, and we have also 
imposed this condition when defining centrals in the Y07 cat¬ 
alog. A number of studies indicate that this is not a good as¬ 


sumption for a substantial fraction of groups (see, e.g., Skibba 
|et al.|[2011[ [Masaki et al.|[2013a[ [Hikage et al.|[^13[ ). Since 
galactic conformity manifests as a similarity of galaxy colours 
in groups, any errors in classifying galaxies as centrals and 
satellites, provided they have a minimal impact on overall 
group membership, will tend to make these populations simi¬ 
lar to each other and are therefore likely to induce conformity¬ 
like effects ( [Campbell et al.|2015[ ). Consequently, it is possible 
that our analysis in Section |5.3| somewhat overestimates the 
strength of conformity (i.e., the value of p). However, the re¬ 
sults of Section [5?^ suggest that our conclusions regarding the 
relative strengths of the conformity signal at various scales are 
qualitatively robust to such errors. It will nevertheless be in¬ 
teresting to revisit this issue in future work using more refined 
criteria for identifying centrals. 


A proper comparison with data will require accounting 
for projection effects both in the red fraction measurements 
as well as the definition of isolated primaries themselves. Al¬ 
though these will dilute the signal, its strength should also in¬ 
crease upon including galaxies with smaller stellar masses (and 
stronger assembly bias) than we could access in our mocks, and 
overall we expect that the difference between large scale 1-halo 
and genuine 2-halo effects will remain measurable. Addition¬ 
ally, we expect that simultaneous measurements of 1-halo and 


large scale conformity together with colour-dependent clus¬ 
tering will be required to break the weak degeneracy between 
halo mass dependence of the central red fraction and the level 
of conformity. Such a joint analysis will also be important from 
the point of view of obtaining unbiased HOD fits that account 
for galactic conformity (Zentner et al.|2014l. 


7 SUMMARY AND CONCLUSIONS 


We have introduced a flexible model of galactic conformity 
within the Halo Occupation Distribution (HOD) framework by 
modifying and extending the algorithm described by |Skibba| 
& Sheth (20091. By construction, our mock galaxy cata¬ 


logs show good agreement with measurements of conformity- 
independent variables in the Yang et al. ( 2007| Y07) group 
catalog based on DR7 of the SDSS. These variables include 
the satellite fraction, the red fractions of all galaxies, centrals 
and satellites, and the satellite quenching efficiency Ssat (Fig¬ 
ure]^, as well as the all-galaxy luminosity and stellar mass 
functions (Figure]^. Galaxy luminosities are assigned using 
the HOD (we use the calibration by Zehavi et al. 20111, colours 


are assigned using colour-luminosity fits to SDSS data (Fig¬ 
ure and equations and stellar masses are assigned using 
a colour-dependent mass-to-light ratio (equation]^, also fit 
to SDSS data. Our mock catalogs are luminosity-complete for 
Mr < —19.0 and mass-complete for logj^Q(m*) > 9.9 (Figuresj^ 
and|^. 

Our algorithm introduces conformity between the colours 
of the central and satellites of a group (a 1-halo effect) by us¬ 
ing a tunable group quenching efficiency p to correlate these 
colours with the concentration of the parent dark halo of the 
group. Halo concentration (which has a scatter at fixed halo 
mass) is therefore identified as the “hidden variable” in our 
model which causes galactic conformity even in halos of fixed 
mass. Halo assembly bias then leads to a 2-halo effect at very 
large scales (this is our default model), which we can also 
switch off by randomizing halo concentrations among halos 
of fixed mass (we call this the no-2h model). The latter is 
a useful toy example in which conformity arises due to some 
unspecified property of halos (e.g., this might be a coupling 
between star formation activity and the hot gas content in a 
halo) that does not exhibit large scale environmental effects. 
We have performed various tests to study the nature of the 
signal, including changing p (Figure in the presence or ab¬ 
sence of 2-halo conformity (Figure [^ for different choices of 
isolation criteria - we used isolated primaries (Figures|^and[^ 
as well as group centrals (Figure |A2| ). Our main results can 
be summarized as follows. 


• We find that setting p = 0.65 gives a 1-halo conformity 
signal in the mocks which agrees well with the corresponding 
signal in the Y07 catalog (Figure]^. The signal manifests as a 
difference between the red fractions of satellites in groups with 
red and blue centrals. Additionally, in the Y07 catalog, the red 
fraction of satellites with blue centrals is remarkably similar to 
the red fraction of centrals at all masses; our mocks correctly 
reproduce this trend down to their completeness limit. 

• The above value of p also leads to a signal at large scales 
in the mocks (Figure]^; specifically, we see a significant dif¬ 
ference between the galaxy red fractions surrounding red and 
blue isolated primaries out to separations < 4Mpc. Interest¬ 
ingly, we find that this signal is dominated by the 1-halo con¬ 
tribution of large halos with logj^Q(m) > 13.25 (Figure |^, 
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and persists even when we switch off 2-halo conformity in our 
no-2h model. This implies that the observation of such a dif¬ 
ference (e.g. Kauffmann et al.|2013 1 is not conclusive evidence 
that galactic conformity arises from halo assembly bias, since 
it could also arise from 1-halo effects “leaking” to large scales 
due to averaging over halo mass (section]^. 

• At even larger scales (> 8Mpc in Figure]^, the signals 
with and without 2-halo conformity do become distinct, with 
the genuine 2-halo signal remaining significantly elevated com¬ 
pared to the global average red fraction out to separations in 
excess of ISMpc (at least in the 3-d case that we consider). 
This 2-halo signal is absent at large halo masses and only 
emerges at smaller masses 11.25 < logiQ(m) < 12.25 (Fig¬ 
ure , being qualitatively consistent with expectations from 
halo assembly bias (section and Appendix |A2[ | . We there¬ 
fore suggest that this elevation compared to the global average 
could be a more robust indicator of large scale galaxy assem¬ 
bly bias than is the difference between red fractions at scales 
< 4Mpc. 

We end with a brief discussion of future extensions of 
our work. In a forthcoming paper (Kovac et ah, in prepara¬ 
tion), we will compare the large scale signal in our mocks, after 
accounting for projection effects, with a corresponding mea¬ 
surement in the Y07 catalog to determine whether or not the 
observed large scale conformity is due to galaxy assembly bias 
or a residual of some other 1-halo process. Tests at fixed group 
richness will also be useful in answering this question. It will 
also be interesting to explore the use of galaxy lensing and/or 
traditional correlation function analyses to validate the con¬ 
nection between galaxy colours and host halo concentrations 
assumed in this work, particularly to check for consistency 
with the value of p presented above. In another publication 
(Pahwa et ah, in preparation), we will present the analytical 
formalism for including conformity in the HOD framework; we 
will show that this requires straightforward modifications in 
existing HOD pipelines. 

Our algorithm can also be extended to include radial pro¬ 


files for satellite velocity dispersions and colours (Prescott 
|et al.|[2011[ [Hartley et al.|[2015p , as well as concentration- 
dependent satellite abundances (Wechsler et al. 2006 Mao 


et al. 20151. The signal at high redshift is interesting too; 


this could be modelled using analytical prescriptions tuned to 
match high-redshift luminosity function and clustering data 


(Tacchella et al. 2013 Behroozi et al. 2013a Jose et al. 
|2014| ). Finally, it will be extremely interesting to extend our 
algorithm to lower stellar masses using accurate faint-end 
(Afr,max ~ —16.0) HOD fits and higher resolution Y-body 
simulations. We estimate that a simulation with 1024^ parti¬ 
cles in a (100/i“^Mpc)® box will allow us to create catalogs 
that are mass-complete for logjg(m*) > 8.7; comparing these 
to data will provide us with stringent tests on the nature of 
the conformity signal. 


ACKNOWLEDGEMENTS 

We thank J. Woo and M. Shirazi for help with the Kcorrect 
tools, and an anonymous referee for an insightful report that 
has helped improve the presentation. We are grateful to the 
SDSS collaboration for publicly releasing their data set, and 
Yang et al. for releasing their updated group catalog. We thank 
O. Hahn, V. Springel and P. Behroozi for making their codes 
publicly available. We gratefully acknowledge computing fa¬ 


cilities at ETH, Zurich and lUCAA, Pune. Our mock catalogs 
are available upon request. 


REFERENCES 

Abazajian K. N. et ah, 2009, ApJS, 182, 543 
Balogh M. L., Morris S. L., 2000, MNRAS, 318, 703 
Baugh C. M., 2006, Reports on Progress in Physics, 69, 3101 
Behroozi P. S., Wechsler R. H., Conroy C., 2013a, ApJ, 770, 
57 

Behroozi P. S., Wechsler R. H., Wu H.-Y., 2013b, ApJ, 762, 
109 

Berlind A. A. et al., 2006, ApJS, 167, 1 
Berlind A. A., Weinberg D. H., 2002, ApJ, 575, 587 
Birrer S., Lilly S., Amara A., Paranjape A., Refregier A., 2014, 
ApJ, 793, 12 

Blanton M. R. et ah, 2003, ApJ, 592, 819 
Blanton M. R., Roweis S., 2007, AJ, 133, 734 
Blanton M. R. et ah, 2005, AJ, 129, 2562 
Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., 
Kravtsov A. V., Klypin A. A., Primack J. R., Dekel A., 
2001, MNRAS, 321, 559 

Campbell D., van den Bosch F. C., Hearin A., Padmanab- 
han N., Berlind A., Mo H. J., Tinker J., Yang X., 2015, 
MNRAS, 452, 444 

Carollo C. M. et ah, 2014, ArXiv e-prints: 1402.1172 
Carretero J., Castander F. J., Gaztanaga E., Crocce M., Fos- 
alba P., 2015, MNRAS, 447, 646 
Cen R., 2014, ApJ, 781, 38 
Cooray A., Sheth R., 2002, Phys. Rep., 372, 1 
Coupon J. et ah, 2015, MNRAS, 449, 1352 
Croton D. J. et ah, 2006, MNRAS, 365, 11 
Dalai N., White M., Bond J. R., Shirokov A., 2008, ApJ, 687, 
12 

Deason A. J., Conroy C., Wetzel A. R., Tinker J. L., 2013, 
ApJ, 777, 154 

Desjacques V., 2008, MNRAS, 388, 638 
Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 
Fabian A. C., Sanders J. S., Allen S. W., Crawford C. S., 
Iwasawa K., Johnstone R. M., Schmidt R. W., Taylor 
G. B., 2003, MNRAS, 344, L43 
Fakhouri O., Ma C.-R, 2010, MNRAS, 401, 2245 
Gao L., Springel V., White S. D. M., 2005, MNRAS, 363, L66 
Gunn J. E., Gott, III J. R., 1972, ApJ, 176, 1 
Hahn O., Abel T., 2011, MNRAS, 415, 2101 
Hahn O., Porciani C., Dekel A., Garollo C. M., 2009, MNRAS, 
398, 1742 

Hansen S. M., Sheldon E. S., Wechsler R. H., Koester B. P., 
2009, ApJ, 699, 1333 

Hartley W. G., Conselice C. J., Mortlock A., Foucaud S., 
Simpson C., 2015, MNRAS, 451, 1613 
Hearin A. P., Watson D. F., 2013, MNRAS, 435, 1313 
Hearin A. P., Watson D. F., Becker M. R., Reyes R., Berlind 
A. A., Zentner A. R., 2014a, MNRAS, 444, 729 
Hearin A. P., Watson D. F., Becker M. R., Reyes R., Berlind 
A. A., Zentner A. R., 2014b, MNRAS, 444, 729 
Hearin A. R, Watson D. F., Bosch F. C. v. d., 2015, MNRAS, 
452, 1958 

Hikage G., Mandelbaum R., Takada M., Spergel D. N., 2013, 
MNRAS, 435, 2345 

Hinshaw G. et al., 2009, ApJS, 180, 225 


© 2015 RAS, MNRAS 000, 






























16 A. Paranjape, et al. 


Hirschmann M., De Lucia G., Wilman D., Weinmann S., 
lovino A., Cucciati O., Zibetti S., Villalobos A., 2014, 
MNRAS, 444, 2938 

Hopkins P. F. et al., 2010, ApJ, 724, 915 
Jing Y. P., Suto Y., Mo H. J., 2007, ApJ, 657, 664 
Jose C., Srianand R., Subramanian K., 2014, MNRAS, 443, 
3341 

Kauffmann G. et al., 2003, MNRAS, 341, 33 
Kauffmann G., Li C., Zhang W., Weinmann S., 2013, MNRAS, 
430, 1447 

Knobel G., Lilly S. J., Woo J., Kovac K., 2014, ArXiv e-prints: 
1408.2553 

Koester B. P. et al., 2007, ApJ, 660, 221 
Kovac K. et al., 2014, MNRAS, 438, 717 
Krause E., Hirata C. M., Martin C., Neill J. D., Wyder T. K., 
2013, MNRAS, 428, 2548 

Larson R. B., Tinsley B. M., Galdwell G. N., 1980, ApJ, 237, 
692 

Lu Y., Yang X., Shen S., 2015, ApJ, 804, 55 
Ludlow A. D., Navarro J. F., Angulo R. E., Boylan-Kolchin 
M., Springel V., Frenk C., White S. D. M., 2014, MNRAS, 
441, 378 

Mao Y.-Y., Williamson M., Wechsler R. H., 2015, ApJ, 810, 
21 

Masaki S., Hikage G., Takada M., Spergel D. N., Sugiyama N., 
2013a, MNRAS, 433, 3506 

Masaki S., Lin Y.-T., Yoshida N., 2013b, MNRAS, 436, 2286 
Moore B., Katz N., Lake G., Dressier A., Oemler A., 1996, 
Nature, 379, 613 

More S., Miyatake H., Mandelbaum R., Takada M., Spergel 
D. N., Brownstein J. R., Schneider D. P., 2015, ApJ, 806, 
2 

More S., van den Bosch F. C., Gacciato M., Skibba R., Mo 
H. J., Yang X., 2011, MNRAS, 410, 210 
Muldrew S. 1. et al., 2012, MNRAS, 419, 2670 
Navarro J. F., Frenk G. S., White S. D. M., 1996, ApJ, 462, 
563 

Navarro J. F., Frenk G. S., White S. D. M., 1997, ApJ, 490, 
493 

Oppenheimer B. D., Dave R., Keres D., Fardal M., Katz N., 
Kollmeier J. A., Weinberg D. H., 2010, MNRAS, 406, 
2325 

Osmond J. P. F., Ponman T. J., 2004, MNRAS, 350, 1511 
Paranjape A., 2014, Phys. Rev. D, 90, 023520 
Peng Y.-j. et al., 2010, ApJ, 721, 193 

Peng Y.-j., Lilly S. J., Renzini A., Garollo M., 2012, ApJ, 757, 
4 

Phillips J. I., Wheeler G., Boylan-Kolchin M., Bullock J. S., 
Gooper M. G., Tollerud E. J., 2014, MNRAS, 437, 1930 
Prescott M. et al., 2011, MNRAS, 417, 1374 
Rawlings S., Jarvis M. J., 2004, MNRAS, 355, L9 
Ross A. J., Brunner R. J., 2009, MNRAS, 399, 878 
Scoccimarro R., 1998, MNRAS, 299, 1097 
Sheth R. K., Tormen G., 2004, MNRAS, 350, 1385 
Skibba R., Sheth R. K., Connolly A. J., Scranton R., 2006, 
MNRAS, 369, 68 

Skibba R. A., 2009, MNRAS, 392, 1467 
Skibba R. A., Sheth R. K., 2009, MNRAS, 392, 1080 
Skibba R. A., Sheth R. K., Croton D. J., Muldrew S. L, Abbas 
U., Pearce F. R., Shattow G. M., 2013, MNRAS, 429, 458 
Skibba R. A., Sheth R. K., Martino M. C., 2007, MNRAS, 
382, 1940 

Skibba R. A. et al., 2014, ApJ, 784, 128 


Skibba R. A., van den Bosch F. C., Yang X., More S., Mo H., 
Fontanot F., 2011, MNRAS, 410, 417 
Somerville R. S., Dave R., 2014, ArXiv e-prints: 1412.2712 
Springel V., 2005, MNRAS, 364, 1105 
Tacchella S., Trent! M., Garollo C. M., 2013, ApJ, 768, L37 
Tinker J. L., George M. R., Leauthaud A., Bundy K., 
Finoguenov A., Massey R., Rhodes J., Wechsler R. H., 

2012, ApJ, 755, L5 

Toomre A., Toomre J., 1972, ApJ, 178, 623 
van den Bosch F. C., Aquino D., Yang X., Mo H. J., Pasquali 
A., McIntosh D. H., Weinmann S. M., Kang X., 2008, 
MNRAS, 387, 79 

van den Bosch F. C., More S., Gacciato M., Mo H., Yang X., 

2013, MNRAS, 430, 725 

van der Wei A., Bell E. F., Holden B. P., Skibba R. A., Rix 
H.-W., 2010, ApJ, 714, 1779 

Wang L., Weinmann S. M., De Lucia G., Yang X., 2013, MN¬ 
RAS, 433, 515 

Wang W., White S. D. M., 2012, MNRAS, 424, 2574 
Watson D. F. et al., 2015, MNRAS, 446, 651 
Wechsler R. H., Bullock J. S., Primack J. R., Kravtsov A. V., 
Dekel A., 2002, ApJ, 568, 52 

Wechsler R. H., Zentner A. R., Bullock J. S., Kravtsov A. V., 
Allgood B., 2006, ApJ, 652, 71 
Weinmann S. M., van den Bosch F. C., Yang X., Mo H. J., 
2006, MNRAS, 366, 2 

Wetzel A. R., Tinker J. L., Conroy C., 2012, MNRAS, 424, 
232 

Wetzel A. R., Tinker J. L., Conroy C., van den Bosch F. C., 
2013, MNRAS, 432, 336 

Wetzel A. R., White M., 2010, MNRAS, 403, 1072 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Yang X., Mo H. J., van den Bosch F. C., Jing Y. P., 2005, 
MNRAS, 356, 1293 

Yang X., Mo H. J., van den Bosch F. C., Pasquali A., Li C., 
Barden M., 2007, ApJ, 671, 153 
York D. G. et al., 2000, AJ, 120, 1579 
Zehavi 1. et al., 2011, ApJ, 736, 59 
Zehavi 1. et al., 2005, ApJ, 630, 1 

Zentner A. R., Hearin A. P., van den Bosch F. C., 2014, MN¬ 
RAS, 443, 3044 

Zhuravleva 1. et al., 2014, Nature, 515, 85 


APPENDIX A: 

Al The (in)dependence of galaxy red fractions on 
halo mass 

As mentioned in the main text, one might expect in general 
that the distribution of galaxy colours depends on both galaxy 
luminosity and host halo mass. We explore this idea further 
in this Appendix. 

Let the all-galaxy red fraction be some function 
p(red|Mr, m) at fixed galaxy luminosity and parent halo mass. 
Separating the red fractions of centrals and satellites, we can 
write this as 

p{xeA\Mr, m) = p(red|cen, Mr, m)p(cen|Mr, m) 

J- p(red|sat, Mr, m)p(sat|Mr, m) , (Al) 

where p(sat|Mr,m) gives the fraction of galaxies that are 
satellites [so that p(cen|Mr,m) = 1 — p(sat|Mr,m)], and 
p(red|sat/cen, Mr, m) gives the fraction of satellites/centrals 
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that are red. Any mass dependence in p(red|sat/cen, Mr, m) 
could lead to “conformity-like” effects in analyses that average 
over halo mass (as we saw in Fignre|^. What is usually re¬ 
ferred to as conformity, however, is the statement that central 
and satellite colours are correlated even at fixed halo mass. 
Since we are trying to model the latter effect, it is important 
to understand the former as well. 

For ease of notation, we define the following quantities 
that are fixed by the HOD: 


A/’sat(< Mr\m) — /cen (< Mr\m)Nsat{< Mr\m) 
A/’aat(Mr|m) = 9A/’aat(< Mr\m)/dMr 
fcen{Mr\m) = 9/cen(< Mr\m)/dMr , (A2) 


where /cen(< Mr\m) and fVaat(< Mr\m) were introduced in 
section pHj The satellite fraction is then also fixed by the HOD 
as: 


p(sat|Mr, m) 


_A/'aat(MT.|m)_ 

/cen (Mn|m) -I-A/’aat(Mr|m) ’ 


(A3) 


In principle, however, we must calibrate the luminos¬ 
ity and mass-dependence of both p(red|sat, Mr, m) and 
p(red|cen, Mr, m). The S09 algorithm simplifies the problem 
by postulating that these two functions are approximately in¬ 
dependent of parent halo mass. In other words, the S09 model 
states 


p(red|Mr, m) = p(red|cen, Mr) p(cen|Mr, m) 

-I- p(red|sat, Mr) p(sat|Mr, m), 


(A4) 


so that the all-galaxy red fraction still inherits a mass de¬ 
pendence from the HOD. S09 proceed to calibrate the func¬ 
tion p(red|sat, Mr) using clustering measurements of SDSS 
galaxies, and fix p(red|cen, Mr) by demanding that the halo- 


averaged version of equation (A4| return the measured all¬ 


galaxy red fraction as a function of luminosity: 

p(red|Mr) = p(red|cen, Mr)p(cen|Mr) 

-I- p(red|sat, Mr) p(sat|Mr 


where 


p(sat|Mr) = 


/ dmn(m) A/'sat(Mr|m) 


/ dmn{m) [/cen(Mr|m) -|- A/’sat(Mr|m)] 


(AS) 


(A6) 


and p(cen|Mr) = 1 — p(sat|Mr), with n{m)dm being the num¬ 
ber density of m-halos. 


The assumptions underlying equation (A41 were tested by 
|Skibba| ( |2009| ) by comparing colour distributions at fixed group 
richness as predicted by the S09 model with those in the Y07 


and Berlind et al. (20061 catalogs. Overall the data appear con¬ 


sistent with the assumption that the red fractions at fixed Mr 
- of both satellites and centrals - do not depend on halo mass. 
As S09 point out, this may be expected for satellites, whose lu¬ 
minosity function is observed to be approximately halo mass- 
independent ( Skibba et al.|20OT Hansen et al.|[2009 l; the sit¬ 
uation for the centrals, however, is less clear since their lumi¬ 
nosity function does depend strongly on halo mass, at least 
for small masses. Also, the calibration by S09 using clustering 
measurements ignored effects of conformity and is therefore 
difficult to interpret in the present context. It is conceivable 
that a weak halo mass dependence p(red|cen, Mr, m) might 
have litte consequence for the fixed-richness tests of [Skibbli] 
( |2009[ ) but might have noticeable effects, e.g., in the red frac¬ 
tion profiles of Figures and 

To assess the importance of such effects, we explore an 
alternative toy model which also assumes a mass-independent 



logio(m, / h"^Mo) 


Figure Al. Comparison of the conformity-independent observ¬ 
ables /red|cen (green circles) and Ssat (magenta squares), and the 
conformity-dependent /red|sat,rc (red triangles) and /red|sat,bc (blue 
inverted triangles) in the Y07 catalog (open symbols with dashed 
lines) and in mocks that used p = 0.65 and a halo mass dependent 
central red fraction given by equations -| |A9[ | (filled symbols 
with solid lines). The error bars on the mock results show the r.m.s. 
fluctuations around the mean value over 10 realisations. 


satellite red fraction, but has a mass-dependent central 
red fraction p(red|cen, Mr, m) whose average satisfies equa¬ 
tion (A5l, with p(red|sat, Mr) and p(red|Mr) given by equa¬ 
tion (121 and the first of equations respectively. At fixed 
luminosity, this model preferentially places red centrals in 
massive halos. In particular, denoting prc = p(red|cen, Mr), 
we use 


p(red|cen, Mr, m) = pr, 


g{Mr, m) 
(p(Mr,m)) 


where we set 


p(Mr, m) = Prc + (1 - Prc) 


tanh(m/Mniin) 

tanh(l) 


(A7) 


(AS) 


with 


J dmn(m)/cen(Mr|m) 


and where Mmin(Mr) is the mass scale where /cen(< Mr|m) = 
1/2. The model is chosen such that in the limit of a step¬ 
like function fcen{< Mr|m) = 0(m — Mmin(Mr)), we would 
have (gr(Mr,m)) = 1 and p(red|cen, Mr, m = Mmin) ~ 
p(red|cen, Mr). The HOD we use has a smoother mass de¬ 
pendence in /cen(< Mr|m) and therefore our model requires 
a renormalisation by {g{Mr,m)}. The model ensures 0 < 
p(red|cen, Mr, m) < 1 in practice for all Mr and m. 

Figure ED shows the conformity-dependent satellite red 
fractions as well as the central red fraction and Ssat when set¬ 
ting p = 0.65 in this model as the filled symbols joined by 
solid lines, formatted as in Figures]^ andThe open symbols 
with dashed lines again show corresponding measurements in 
the Y07-Mass catalog, and we see a good agreement between 
the mocks and the data. The 2-halo signal predicted by this 
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Figure A2. Same as the right panels of Figures and using the same set of mocks, except that the galaxy red fractions are measured 
centered on group centrals (centrals with at least one associated satellite) rather than isolated primaries. See text for a discussion. 


model is shown in Figure!^ as the open symbols/dashed lines 
and is discussed in section[6l 


A2 Large scale trends with and without conformity 


Here we further explore the nature of large scale conformity in 
our mocks and argue that the trends seen in Figures and 
can be understood as a consequence of the relative bias of blue 
and red galaxies, together with the nature of halo assembly 
bias. 

It is interesting to compare the results in the right pan¬ 
els of Figures and with similar measurements performed 
around “group centrals” rather than isolated primaries, where 
a group central is a central with at least one associated satel¬ 
lite. Firstly, we find that the number density of group cen¬ 
trals in our simulation box (~ 1.1 x 10“^(/i/Mpc)®) is only 
about 16% that of isolated primaries, meaning that isolated 
primaries predominantly comprise of isolated singletons which 
occupy low mass halos, which is not surprising considering that 
~ 90% of all centrals are singletons. The all-galaxy red fraction 
surrounding group centrals is shown in Figure [A^ whose left 
and right panels are formatted exactly like the right panels of 
Figures and respectively. We see, e.g., that the magnitude 
of the difference between the p = 0.99 red fractions around 
red and blue group centrals is substantially reduced as com¬ 
pared to that around isolated primaries. At scales larger than 
~ 15Mpc, the no-2h mocks show red fractions consistent with 
the global average, while the default mocks with increasingly 
larger p show correspondingly larger red fractions. 

To understand these trends, consider the simpler mea¬ 
surement of the galaxy red fraction at a distance r from an 
arbitrary red central, /red|rc(>")- A. straightforward calculation 
shows that we can write this as 


/red|rc(r) 


f. I / nb(l-f^b,rc(r-)) \\ 


(AlO) 


where hb(r) is the mean number density of blue (red) galaxies 
and ^b(r),rc(r) is the cross-correlation function between blue 
(red) galaxies and red centrals. A similar expression holds for 


the red fraction around blue centrals /red|bc(?'), with the re¬ 
placement ‘rc’—^-'bc’ in the correlation functions on the r.h.s. 
of equation ( |A10[ ). 

The correlation functions above can be split into contri¬ 
butions arising from galaxies in the same halo as the central 
object or in different halos (the 1-halo and 2-halo terms, re¬ 
spectively); ^(r) = (r)(r). Recall that 1-halo terms 

dominate at small separations and 2-halo terms at large sepa¬ 
rations. For sufficiently large separations the 2-halo terms also 
become small compared to unity. We now observe the follow¬ 
ing: 


• Red fraction profile: 

At scales large enough that the 2-halo contributions dominate 
but are not yet negligible compared to unity, /red|rc(f') will 
have an r-dependence given by replacing ^(r) —>■ in 

equation (AlO I. In general, the correlation functions would not 
conspire to precisely cancel the r-dependence in this regime. 
Therefore, only when both 2-halo correlations have become 
substantially smaller than unity will the asymptotic value 
/red|rc(r) —^ (1 -|- fih/hr)~^ = /red be reached. To get a rough 
estimate of when this will occur, let us ignore complications 
introduced by mass completeness, colour cuts, etc. and simply 
use a power law fit for the all-galaxy correlation functio n of 
our luminosity complete sample: ^(r) ~ (r/7Mpc)“^'®® (Ze- 


havi et al. 20111. It is therefore not surprising that all our 


data sets, including the ones with p = 0.01, show significant 
red fraction profiles at distances less than ~ 15Mpc. 

• /red|rc(r) = /red|bc(r) at large r for all models: 

At large separations, it is reasonable to assume that halo 
bias is approximately linear and scale-independent. The 2-halo 
cross-correlations above can then be written as Fourier trans¬ 
forms of cross-power spectra that take the form RbW^rc(^) ~ 

Riin(fe) &rc feb(r)(fc) for CbM',rc(’')> "here 6b(r)(fc) is the number- 
weighted halo bias associated with all blue (red) galaxies while 
6 rc is the corresponding term for red centrals alone, and Piin(fc) 
is the dark matter power spectrum in linear theory (see, e.g., 


Cooray & Sheth 20021. (The scale dependence in the all-galaxy 


bias arises from the profiles of satellites within a halo, while 
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distance from isolated primary (Mpc) 


Figure A3. Same as Figure]^ showing the large scale conformity signal for the default and no-2h mocks around isolated primaries with 
9.9 < logjQ(m*) < 10.1 (left panel) and logjQ(m*) > 10.4 (right panel). 


the central-only term has no such scale dependence.) In the 
regime where both these terms are larger than unity, the scale 
independent 6rc cancels in their ratio. A similar factor febc can¬ 
cels when constructing /red|bc(t') as well, which means that we 
will have /red|rc(f) ~ /red|bc(f) iu this regime for any value of 
p. Indeed, we see this behaviour at large scales in Figure [A^ 

• All red fractions are larger than /red: 

The red fractions (even in the absence of 2-halo conformity) 
approach the global mean from above rather than below be¬ 
cause the average bias factors of blue and red galaxies satisfy 
6 b < 6r. This may seem counter-intuitive since, on average, 
blue galaxies in a mass-complete sample are more luminous 
and therefore preferentially occupy more massive halos than 
red galaxies. However, the halo mass distribution of red galax¬ 
ies is substantially broader than that of blue galaxies, with a 
longer tail at high masses that pushes the number-weighted 
red galaxy bias 6r to be larger than bb. In the luminosity com¬ 
plete sample, this effect is further enhanced due to the inclu¬ 
sion of a large number of faint blue objects (predominantly 
centrals) which also live in low mass halos. 

• Increase of 2-halo conformity red fractions with p: 

The reason that the red fractions in the presence of 2-halo con¬ 
formity increase with p is the presence of halo assembly bias 
at large scales. In this case the bias factors above contain inte¬ 
grals of the form ~ f dlnm f dsp(s|m) 6(m, s) ‘I>biue(red) (ui-, s) 
which couple the concentration-dependence of halo bias to 
that of the red fraction in our model. The number weighting 
of bias means that small mass halos (m 2 oob ~ ~ 

0.25M,) will dominate the integral. If the bias b(m, s) for such 
halos were monotonically increasing with s, then a straightfor¬ 
ward calculation shows that a large positive value of p would 
enhance the difference between 6b and 6^ discussed previously, 
consequently increasing the red fractions. Although the actual 
behaviour of the bias for masses ~ 0.25M* is non-monotonic 
and shows a minimum around s = 0 (see, e.g.. Figure 4 of 

, the fact that the majority of galaxies 
sample are red for any halo mass means 
that this trend does not change. Interestingly, the amplitude 


of this signal appears stable against a possible weak mass de¬ 
pendence in the central red fraction (see Figure]^, suggesting 
that an enhanced value of /red|rc(r) — /red|bc(r) at r > 8Mpc 
as compared to the global average could be a robust indicator 
of 2-halo conformity. 

• /red|rc(7') > /red|bc(»') for v < 5Mpc for p = 0.99 with 
2-halo conformity: 

For p = 0.99, the term “red fraction” is essentially the same as 
“fraction with large halo concentration s” (equation]^. When 
computing red fractions we are then actually asking: what is 
the fraction of large-s objects around (a) large-s objects and 
(b) small-s objects? If we now focus on separations r < 5Mpc, 
we are directly studying a version of halo assembly bias in 
a range of scales (the so-called “l-halo to 2-halo transition” 
regime) where the previous assumptions regarding linear scale- 
independent biasing are no longer valid. 

Consider centrals in low mass halos first (e.g., the isolated 
singletons from earlier). Due to assembly bias, the ones with 
large s will cluster strongly with other large-s centrals of sim¬ 
ilar halo mass, while the small-s centrals will be close to un¬ 
biased. The number weighting will ensure that this behaviour 
gives the dominant effect in the vicinity of small halo centrals, 
with little contribution from satellites and centrals in more 
massive halos. This will lead to a large difference in /red|rc(’') 
and /red|bc(’') ut Small distances, which must then smoothly 
match the large scale trend of /red|rc(r) ~ /red|bc(r'), consistent 
with what is seen in Figures]^ and 

Group centrals that reside in more massive halos are likely 
to have other high mass halos nearby due to larger average 
halo bias, so that assembly bias will still be active. However, 
contamination due to nearby singletons (which reside in less 
clustered, lower mass halos) and the fact that assembly bias 
reverses its trend around the global characteristic halo mass 
will mean that the strength of the difference in red fractions 
will be smaller here. This is consistent with the results in Fig¬ 
ure El 

• p = 0.99 without 2-halo conformity: 

For the model with p = 0.99 and 2-halo conformity switched 


Wechsler et al. 2006 


in the mass-complete 
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11 12 13 14 15 

log(mhaio/ h'Mo) 


Figure A4. Joint distribution of stellar mass of centrals and their 
parent halo mass. For comparison, the smooth yellow curves show 
the median stellar mass at fixed halo mass (solid curve) and its 
1-sigma scatter (dashed curves) as calibrated by jBehroozi et al.| 
| |2013a[ | at redshift 2 : = 0. The horizontal dotted line shows our 
mass completeness limit at logjQ(m*) = 9.9. 


off, we expect that at large scales /redpc = /red|bc while at 
small scales /red | re 1 and /red | be 0. The transition be¬ 
tween these extremes will be dramatic but smooth due to av¬ 
eraging over halo mass, and should occur in the same “l-halo 
to 2-halo transition” regime discussed above. This is consistent 
with what we see in the right panel of Figurej^ where the tran¬ 
sition occurs around r ~ 3-4Mpc. In the left panel of Figure [X^ 
we essentially see only the large scale behaviour; presumably 
this is because we do not impose any isolation criterion in 
the group central measurement and therefore end up selecting 
very small groups (of which there are many). The isolation 
criterion in Figure]^ on the other hand would select relatively 
large groups (in addition to the very high number of single- 
tons which would however not contribute any 1-halo signal). 
This effect adds to and confuses the transition regime 2-halo 
conformity effect discussed above, particularly when studying 
projected signals. We discuss this further in section 

Finally, the two panels of Figure [A^ are formatted iden¬ 
tically to Figure and show that the large scale conformity 
signal in our mocks is essentially independent of the stellar 
mass of the isolated primary, an effect we attribute to the 
large scatter in the distribution of halo masses at fixed stellar 
mass as seen in Figure [Xd] 
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